Observation of anisotropic effect of antiferromagnetic ordering on the 
superconducting gap in ErNi2B2C 



N. L. Bobrov *, V. N. Chernobay, Yu. G. Naidyuk, L. V. Tyutrina, I. K. Yanson 
B.I. Verkm Institute for Low Temperature Physics and Engineering, 
NAS of Ukraine, J^l, Lenin Prospect, 61103, Kharkiv, Ukraine 

D. G. Naugle and K. D. D. Rathnayaka 
Department of Physics Texas A&M University, College Station TX 77843-4242, USA 

(Dated: July 6, 2010) 

The point-contact (PC) spectra of the Andreev reflection dV/dl curves of the superconducting 
rare-earth nickel borocarbide ErNi2B2C (Tc ~11 K) have been analyzed in the "one-gap" and 
"two-gap" approximations using the generalized Blonder-Tinkham-Klapwijk (GBTK) model and 
the Beloborod'ko (BB) model allowing for the pair-breaking effect of magnetic impurities. Exper- 
imental and calculated curves have been compared not only in shape, but in magnitude as well, 
which provide more reliable data for determining the temperature dependence of the energy gap 
(or superconducting order parameter) A(T). The anisotropic effect of antiferromagnetic ordering at 
Tjv «6 K on the superconducting gap/order parameter has been determined: as the temperature 
is lowered, A decreases by ~25% in the c-direction and only by ~4% in the afe-plane. It is found 
that the pair-breaking parameter increases in the vicinity of the magnetic transitions, the increase 
being more pronounced in the c-direction. The efficiency of the models was tested for providing 
A(T) data for ErNi2B2C from Andreev reflection spectra. 

PACS numbers: 74.45.-l-c, 74.50.-|-r, 74.70.Dd 



INTRODUCTION 

Quaternary intermetallic nickel borocarbides (here- 
after borocarbides) of the i?Ni2B2C type (i? is a rare- 
earth element) attract special interest (see surveys pj, 
0, Q and further references) as they include compounds 
with rather high superconducting transition tempera- 
tures (up to Tc ~17 K, i?=Lu) and compounds with dif- 
ferent types of magnetic ordering that include states with 
commensurate and incommensurate spin-density waves. 

Borocarbides have a body centered tetragonal crys- 
talline structure with the ratio c/a~3 P, They have 
a rather complex Fermi surface (FS) consisting of several 
sheets 0, Si • anisotropic f?, '5] and has two char- 

acteristic groups of electrons possessing different Fermi 
velocities vp (this was mentioned in [6| from the de Haas- 
van Alphen experiments 0, [1] ■ The Tc of borocarbides 
is determined not by the total density of states N{Ef), 
but by the contribution to the density of states which 
is made by the slow electrons of the nodal regions 2]. 
In the normal state the transport properties of boro- 
carbides, in particular their resistivity p, are practically 
isotropic [H, ll| because they are related to the groups 
of electrons that have relatively high velocities i^p with 
lower anisotropy and are unrelated to the nodal points 
at the Fermi surface [2]. 

In i?Ni2B2C (i?=Dy, Ho, Er, Tm) compounds the ele- 
ment R contains 4f-electrons with partially filled f-shells 



having a magnetic moment. As a result, Tc of these 
compounds is appreciably lower in comparison with non- 
magnetic borocarbides {R=Y, Lu P,[l,yl). The object of 
this study, ErNi2B2C, undergoes a superconducting tran- 
sition at Tc ~11 K [l|, 0, 0] and two magnetic transitions 
below Tc, which do not destroy superconductivity. The 
AFM ordering occurs at the Neel temperature Tn ~6 K 
when the Er ions form a transverse-polarized incommen- 
surate spin-density wave state [H, 0, 0|- The AFM or- 
dering entails structural distortions and thus reduces the 
crystal symmetry from tetragonal to orthorhombic Q. 
This magneto-elastic effect is regarded as a structural 
Jahn- Teller transition The modulation wave vec- 

tor of the spin-density waves is practically independent 
of temperature and is along the a-axis (direction [100]), 
or the equivalent 6-axis (direction [010]), the spins being 
aligned along the 6-axis or the a-axis, respectively. As 
the temperature lowers further, the compound changes 
into a weakly ferromagnetic state with T^fm ~2.3 K in 
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which a spontaneous vortex lattice is formed 11 1. 

To understand the features of the superconducting 
state in magnetic borocarbides, it is essential to have 
information about the superconducting gap A (magni- 
tude, behavior, anisotropy, etc.) or the superconducting 
order parameter (OP). The investigations of the gap A 
inErNizBzC [H, E [3, El, [H give A = 1.6-^-1.82 meV. 
Its temperature dependence corresponds on the whole to 
the BCS theory It is noted [H, Hj that on a paramag- 
netic - AFM transition A decreases and the pair-breaking 
parameter F [iTl j reaches a maximum in the transition re- 
gion. The authors [l3| interpreted the results using the 
tlieory 0] which predicts a decrease in A due to spin- 
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density wave gaps that open in some parts of the FS. On 
the other hand, the influence of the AFM transition on 
A was not observed in subsequent tunnel measurements 
on ErNiaBjC 



151. 



Detailed point-contact (PC) spectra of the Andreev 
reflection in ErNi2B2C in two principal crystallographic 
directions have been obtained in our recent study |19| . 
The analysis of these spectra shows that: 

1. They are essentially anisotropic and their behavior 
differs qualitatively from that in LuNi2B2C; 

2. The AFM ordering lowered the gap; 

3. The two-gap model can be efficient at describing 
the experimental results. 

However, the large (six) number of fitting parameters in- 
volved in the two-gap model casts some doubt on the 
uniqueness of results. On the other hand, it is quite 
appropriate to determine to what extent the one-gap ap- 
proximation can account for the nontrivial behavior of 
the superconducting OP (gap) In this study we 

have also analyzed the one- and two-gap approximation 
within the Beloborod'ko (BB) model |l9l |20| allowing 
for the pair-breaking influence of magnetic ions and on 
the basis of the generalized Blonder- Tinkham-Klapwijk 
(GBTK) model 2l| allowing for the Dynes pair-breaking 
parameter (l7j . 



EXPERIMENTAL TECHNIQUE 

Here we detail the experimental technique, which was 
described in [3| only briefly because of the space limi- 
tations. The PC measurements were made on ErNi2B2C 
single crystals (Tc «11 K) grown from the melt (Ames 
Laboratory, Prof. P. Canfield's group) and were simi- 
lar to those used in The crystals were thin (0.1 ~ 
0.2 mm) plates with the c-axis normal to the plane of the 
plate. The surface pretreatment was similar to that for 
LuNi2B2C 2^, 24 1 (either by etching in a 5% nitric acid- 
alcohol solution or by of cleavage). The other electrode 
was a high-purity Ag rod or Ag wire of 0.15 mm in diam- 
eter. In the latter case the wire surface was pre- washed in 
concentrated nitric acid. The contact was made between 
the pretreated single crystal surface and the loop-shaped 
wire. The use of the wire loop as a damper improved the 
mechanical stability of the contacts and made it possible 
to measure the characteristics of one of the contacts (in 
the a6-plane) both on heating from 1.45 K to T >11 K 
(normal state) and on subsequent cooling from the nor- 
mal sate to the starting lowest temperature. 

The temperature measurements were made using a 
continuous flow He cryostat (its analog is described in 
[ist). An insert was placed inside the cryostat, which 
made a PC by touching the sample surface with the sil- 
ver electrode at helium temperature. The typical PC 



resistance varied from several to tens of Ohms. For more 
elaborate investigations, the PCs were selected, which 
had the highest possible "tunneling" characteristics seen 
as an intensive maximum in the dV/dI{V) curve at y = 
and the strongest nonlinearity corresponding to at least a 
10% change in dV/dI{V) in the interval ±8 mV. On some 
contacts the dV/dI{V) spectra were measured in the in- 
terval from T,nin=1.45 K to temperatures 1 2 K higher 
than Tc «11 K. The results were quite reproducible irre- 
spective of the resistance of a particular contact. There- 
fore, here we analyze the measurements on two contacts 
along the principal crystallographic directions along the 
c-axis and in the a5-plane. The detailed series of the 
dV/dI{V) spectra were obtained at approximately equal 
temperature intervals. 30 curves were taken at rising 
temperature and 34 curves were measured at lowering 
temperature in the a5-plane, while 49 curves were mea- 
sured in the c-direction at rising temperature. The con- 
tacts remained stable during the whole period of mea- 
surement. 



RESULTS AND DISCUSSION 

The temperature series of dV/dI{V) curves measured 
on two ErNi2B2C— Ag contacts in the a6-plane and along 
the c-axis are shown in Fig.[T] Parameters characterizing 
the contacts and ErNi2B2C are described in [l3|. The 
most important of them is the PC size (diameter) esti- 
mated as dab ~9.1 nm, dc «4.4 nm for the corresponding 
directions (Fig.[T|). The coherence length ^ in this com- 
pound increases from 15 to ~23 nm ^26j when T increases 
from 3 K to 8 K, which satisfies the requirement < ^ 
for the theory [27l |. 

For better visualization, some dV/dI{V) curves of 
Fig. [T] were symmetrized and then scaled by dividing by 
dV/dI{V) in the normal state at the lowest temperature 
and the bias interval ±8 mV. They are shown in Fig. [2] 
along with the results (liires) of the one-gap calculation 
within the GBTK model [2l| . The scaled curved are re- 
duced to equal amplitude ^norm = A{T)/M{T). Here 
Af(T) = A{T)/A{0), is a coefficient, where A{T) is the 
amplitude at the temperature T. A{0) is the amplitude 
of the curve normalized to the normal state at the lowest 
temperature in the bias interval ±8 mV (^norm = ^(0))- 
Note that the smooth jump-free dependence M(T) in the 
a6-plane and c-direction (Fig. [3]) is indicative of the tem- 
perature stability of the point contacts. 

Of interest is the unusual temperature dependence 
of the distance between the minima in dV/dl of the 
ErNi2B2C-Ag point contact (Fig.g]). It differs dras- 
tically from the corresponding dependence in dV/dl of 
a LuNi2B2C—Ag point contact having a similar tunnel- 
ing parameter z [illil. It is known that at low tem- 
peratures the half-distance between the minima in the 
dV/dI{V) curve of high-tunneling [Z ^ \ S-N point con- 
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Voltage, mV 

FIG. 1: Differential resistances of ErNi2B2C-Ag point contact 
formed in two principal directions at different temperatures. 
To avoid overloading, only a few of the curves are shown 

tacts correlates quite well the superconducting energy 
gap [ii]. It should be noted that in the ErNi2B2C-Ag 
PC the two-minima structure of the dV/dl curve persists 
up to Tc in the a6-plane and 0.95 Tc in the c-direction. 
This is typical of tunnel contacts (e.g., [l^l) and rather 
unusual in PCs, where the Andreev reflection is impor- 
tant and the tunnel parameter is Z < 1(0.8). 

Besides, there is another feature in the ErNi2B2C spec- 
tra that is unobservable in LuNi2B2C: the distance be- 
tween the minima in the dV/ dl curve increases with tem- 
perature up to a maximum slightly above the tempera- 
ture of the AFM transition. It is reasonable to attribute 
this behavior to the magnetic transition in ErNi2B2C. 
Such transitions are absent in LuNi2B2C. 

It is also important that the local critical temperature 
in the investigated PCs (at which the main minimum 
disappears from the dV/dl curve) practically coincides 
with Tc of the crystal, which suggests that the properties 
of the material remain unaltered in the contact. 

As is mentioned in the introduction, experimental data 
were analyzed in the one- and two gap approximations 
using two models: 

1. The traditional GBTK modelEU, which includes 
the broadening parameter F [iTf characterizing in- 
elastic pair-breaking processes; 
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FIG. 2: Symmetrized curves of Fig.[T] normalized to the nor- 
mal state at different temperatures. Points show experimental 
data. Lines are theoretical one-gap GBTK calculation. For 
visualization, all the curves are reduced to the same amplitude 
(A(0)=Anorm, See the text) 



2. The BB model [l^l which introduces the pair- 
breaking parameter 7 to account for the finite life- 
time of Cooper pairs due to the pair-breaking ac- 
tion of (disordered) magnetic moments (in our case 
the Er ions possessing a magnetic moment). 

The terms "energy gap" (GBTK model (U) and "order 
parameter" (20| used in this model are of equivalent phys- 
ical sense (see the detailed discussion in [20|, p. 014512-3 
and are therefore denoted identically with A. This, how- 
ever, does not refer to the term "energy gap Aq" in the 
BB model [20|] which differentiates the superconducting 
order parameter A and the energy gap Aq 20]. The 
energy gap Aq and the order parameter A are related as 

Ao-a(i-7^/=^)'^' (1) 

Here 7=l/rsA is the pair-breaking parameter, Ts is the 
electron mean free time under spin-flip scattering. When 
this scattering is absent, Tg tends to infinity and the equa- 
tion describing the current-voltage characteristics (IVCs) 
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FIG. 3: Temperature dependencies of amplitude coefficients 
M{T) obtained from the plots in Fig. [2] and scaling factors 
5 calculated in the one-OP BB approximation and charac- 
terizing the intensity of experimental curves in comparison 
with theoretical one (see Figs. 11, 12) in the afo-plane and the 
c-direction. Here and in the subsequent figures the measure- 
ments in the a6-plane are marked with solid symbols (rising 
temperature, +) and empty symbols (subsequent cooling, - 
). For visualization, a polynomial fit is drawn through the 
points. 



[20| coincides with the corresponding equation of the clas- 
sical BTK theory • The equations describin g; th e IVCs 
of PCs within this model are presented in 2^ 24 [. 
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S, which characterizes the intensity ratio between ex- 
perimental and theoretical curves, must be equal to 1. 
Sometimes this requirement of the GBTK model is vio- 
lated and we have S ^1. Such factors and the ways of 
their selection are considered in Appendix B. 



Calculation in GBTK and BB models with fixed S 
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ONE-GAP APPROXIMATION 

The calculation technique of the most popular GBTK 
model minimizing the r.m.s. deviations F between the 
shapes of experimental curves (see [2Jj, Fig. 3) faces a 
certain problem: when the magnitudes of T and the gap 
A become comparable, the error curve for A has no dis- 
tinct minimum, which is most typical of the c-direction 
(see Appendix, Fig. 18). 

Note that in this comparison of the theoretical and 
experimental curves the parameter F characterizes only 
the degree of their discrepancy in shape, while the dis- 
tinctions in intensity are compensated using a scaling 
factor S = {dV/dI)cxp/{dV/dT)tcor- The scaling factor 



The dependencies A(T) in the a6-plane and c-direction 
calculated with properly chosen S'-factors are shown in 
Fig.m 

Table 1 contains A(0) estimates (bold type) for both 
directions obtained with proper S'-factors. Also, it in- 
cludes GBTK data (A(0)) for the maximum (S=l) and 
minimum (5=0.25) possible cases. Such S values appear 
because in this model the error in F changes only slightly 
when S deviates from its proper value (Figs. [T51B(I)) . 

Note that in the a6-plane the relation 2A/kTc = 3.52 
agrees with the BCS theory at 5=0.31 (our proper selec- 
tion). In the c-direction the BCS relation 2A/kTc — 3.53 
(A(0)) =1.62 meV in the PM region is achieved at 
5=0.54. In this case the error in F increases only slightly. 
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FIG. 5: Temperature dependencies of gaps (GBTK-triangles) 
and OPs (BB-circles) calculated for the best-selected scaling 
factors S (Fig.[20l) for the PC in Fig.[Tl Solid lines are BCS 
extrapolations. 



see Fig.lini)- 



TABLE I: Superconducting energy gaps (GBTK) or OPs 
(BB) A at different scaling factors 5* in the a6-plane and c- 
direction. 



s 


Direction, Face 


GBTK 


BB 


A(0) 


2A/fcTc 


A(0) 


2A/fcTc 


0.25 


ab, PM 


1.8 


3.76 


1.95 


4.08 


c, PM 


2.17 


4.73 


2.36 


5.14 


0.31 


ab, PM 


1.68 


3.52 




0.65 


c, PM 


1.48 


3.23 


1 


ab, PM 


1.03 


2.15 


c, PM 


1.23 


2.68 



The calculation in the BB model gives 2A/kTc = 4-^5. 
This correlates with the tunnel investigation on nonmag- 
netic YNi2B2C, in which 2A/kTc = 5.2 for the maximum 
gap (see Fig.6). 

The critical temperature extrapolated to the paramag- 
netic region of the BCS curve is close to the values of the 
bulk compound Tc==10.64 K (c-direction) and Tc=ll.l K 
(a6-plane ). At the same time at the AFM-PM tran- 
sition the growth of A is essentially dependent on the 



direction. According to the BCS extrapolation, A in- 
creases by 25-28% in the c-direction and only by 4-5% 
in the a6-plane. Proceeding from the magnetic structure 
of ground-state ErNi2B2C [if, this anisotropic influence 
of the magnetic transition on A can be attributed to 
orientationally-dependent spin-density waves. The AFM 
incommensurate ordering below the Neel temperature in- 
duces spin-density waves whose propagation vector q is 
in the a6-plane. Such waves reduce the superconducting 
gap for the electrons having the wave vector k perpen- 
dicular to the vector q, i.e. for the c-direction, due to the 
pair- breaking exchange field diim. The same approach 
was used to interpret the anisotropy of the superconduct- 
ing energy gap in the AFM heavy- fermion URu2Si2 com- 
pound [32| ■ The anisotropic effect of spin-density waves 
is also evident in the behavior of the parameter M (see 
Fig-El, which has an extremum only in the c-direction. 
It is important that M characterizes the "gap minima" 
intensity of the original dV/dl spectra and is unrelated 
to any theoretical model. 

The effect of spin fluctuations and the AFM molecular 
field on the superconducting gap is determined by the 
sum rule [33| . Their competition dictates whether the 
AFM phase will enhance or suppress the pair-breaking 
processes below Tn- The temperature dependence of the 
superconducting gap was calculated within the Chi-Nagi 
model i3j] (Fig. 3 in [33i|). In the paramagnetic region the 
superconducting gap follows the BCS-dependence and in 
the AFM region (below Tn-) its behavior is determined 
by the interaction between the temperature-dependent 
AFM molecular field and the spin-fiuctuation scattering 
of conduction electrons at both magnetic rare-earth ions 
and nonmagnetic impurities. The molecular field opens 
AFM gaps in the some parts of the FS and destroys the 
superconducting gaps in them. Nonmagnetic impurities 
have no effect on the BCS states of a s-wave superconduc- 
tor, but they attenuate the AFM field effect suppressing 
the pairing states of charge-density or spin-density waves. 
The degree of suppression of the superconducting gap is 
dependent on the single crystal perfection - elastic scat- 
tering assists in restoring superconductivity. 

The temperature dependence of the energy gap Aq = 
A (1-72/3)'^/^ (Fig.© can be obtained clearly within 
the BB model [20| as well. The BCS extrapolations 
on changing to gapless superconductivity in the AFM 
(6.6 K) and PM (9.3-9.6 K) regions correlate well in both 
direction with the BCS extrapolations of Tc for the OP5 
in these regions (Fig. [5]). As in the OP case, the gap in- 
creases on changing to the PM state in the c-direction 
and exhibits a monotonic dependence in the a6-plane. 

Note that the results of this study and [l^ (decreasing 
OPs/gaps in the AFM region) correlate with the tem- 
perature dependencies of the coherence length, the pen- 
etration depth and the critical magnetic fields measured 
in single crystalline ErNi2B2C which has features near 
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FIG. 6: The temperature dependencies of the BB-calculated 

. 3/2 



gap Ao = A (^1 



7 
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for the ErNi2B2C-Ag PC in the ab- 
plane and the c-direction. Sohd Unas are BCS extrapolations. 



Tn: an TV-shaped curve of the coherence length and a 
local minimum in the penetration depth ([35|, Fig.l). 
There is also indirect evidence for the anisotropy of A 
which is based on measurements of the anisotropy of 
the upper critical magnetic field Hc2 (it is known that 
Hc2 ^ ^'"^ ^ A'^). A 3-D fourfold modulation of the up- 
per critical field Hc2 was measured in the field H_Lc at 
T=2 K as a function of the direction in the a&-plane (see 
(35} , Fig. 4, insert). Its shape is similar to the anisotropic 
function of the superconducting energy gap in the model 
proposed in [s^. Note a distinct peak at Tn in the de- 
pendence Hc2 (T) in the field along the c-direction (Fig. 2 
in 



35[), which indicates indirectly that the AFM transi- 



tion reduces the superconducting gap. Fig. [7] and Fig. [5] 
illustrate the pair-breaking parameters F and 7 reduced 
to the same dimensionality: 7 is compared with F/A and 
F with 7 • A. The dependencies have features near the 
Neel temperature - a maximum in the c-direction and a 
smeared shoulder in the a&-plane. In the c-direction the 
pair-breaking parameter increases slightly near the tran- 
sition to weak ferromagnetism (~2 K). In the a6-plane 
such increase is hard to identify because of the scatter of 
experimental points. There is a certain correlation be- 
tween the pair-breaking parameter 7 (hence, increasing 
scattering of superconducting electrons) and the tunnel- 
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FIG. 7: Temperature dependencies of broadening parameters 
r and pair-breaking parameters 7 reduced to the same units 
for the ErNi2B2C-Ag PC in the a&- plane and the c-direction. 
Error bars represent 5% deviations A from the minimum at 
curve of r.m.s. deviations of the shape of the theoretical curve 
from experimental one (see Fig. llSp . Since S=const was ob- 
tained by selecting the broadening (pair-breaking) parameter, 
a comparatively slight deviation of the points from the poly- 
nomial fit, the real deviation is much smaller 



ing parameter Z characterizing the potential barrier or 
the scattering at the N — S boundary (see Fig. [8]). Thus, 
the parameter Z may account for both the elastic scat- 
tering intensity and the spin-flip scattering. It is likely 
that the growth of these parameters near the AFM tran- 
sition is more evident in the c-direction because of the 
anisotropic influence of the spin-density waves. On ap- 
proaching Tc , 7 and Z increase in both the afc-plane and 
the c-direction. 



Calculation with a varying S- factor (BB model). 

Note that the use of a fixed S in the BB model re- 
duced the quality of fitting of theoretical to experimental 
curves. The reasons for such reduction in quality are con- 
sidered in Appendix B. The temperature dependencies of 
OPs having freely varying S factors are shown in Fig.[S] 
for the a&-plane and c-direction (two-gap approximation 
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FIG. 8: Temperature dependencies of the pair-breaking parameter 7 and the relative broadening parameter F/A (upper 
pictures) in comparison with the temperature dependencies of the tunnel parameter Z (lower pictures) for the ErNi2B2C-Ag 
PC in the afe-plane and the c-direction. The dependencies 7(r) and F/A(T) are shown on an enlarged scale in the c-direction 
to visualize the maximum near the magnetic transition. The error bars are found as in Figs. [7] 



with free S was used in [19j). In both directions the 
OPs deviate from the BCS curve on approaching Tc (at 
T ~10 K) and turn zero at Tc of the contact. The criti- 
cal temperature obtained for these parameters through a 
BCS extrapolation is Tc=11.9 K in both directions. The 
temperature dependencies of the scaling factors S (Fig. [3]) 
correlate in shape with those of the amplitude coefficients 
M{T). Such temperature dependencies are due to the 
density of states curve in this model which differs from 
the BCS dependence (see the Appendix), as well as to the 
two-gap character of superconductivity in this compound 
whose OPs have significantly different magnitudes in the 
paramagnetic region (detailed in the following section). 

Finally, Fig. [10] illustrates the temperature dependen- 
cies of the pair-breaking parameter 7. On the whole, 
they are similar to the temperature dependencies of the 
pair- breaking parameters for a large OP in the two-gap 
approximation lo'] and differ considerably from the cor- 
responding dependencies in the one-gap approximation 



obtained with a fixed scaling factor (Fig. [5]), which is par- 
ticularly evident in the high-temperature region. This 
may be because in the latter case we try to hold the scal- 
ing factor invariable at the expense of a certain departure 
from coincidence of the shapes of experimental and the- 
oretical curves. In the vicinity of magnetic transitions 
pair-breaking parameters also increase in both a6-plane 
and c-direction. 

Thus, the behavior of OPs and other parameters es- 
timated in the one-gap approximation within the BB 
model with a free factor S is similar qualitatively to the 
results obtained in the two-gap approximation |19| . 



TWO-GAP APPROXIMATION 

The two-gap approximation assumes that the total 
conductivity is a superposition of conductivities from two 
region (bands) of the FS with corresponding gaps. This 
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FIG. 9: Temperature dependencies calculated in the one OP 
approximation (BB model [20j]) with S 7^ const. Solid lines 
are BCS extrapolations 
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FIG. 10: Temperature dependencies of pair-breaking param- 
eters 7 in the one OP approximation with S T^const 



can be expressed for dV/ dl as 



dT ~ d.^ 



s 

If (Ai, 71, Z) + If (A2, 72, ^) (1 - K) 



(2) 



Here the coefficient K accounts for the contribution to 
the conductivity from the FS region with a smaller gap 
Ai, Z is the tunnel parameter, S is the scaling fac- 
tor characterizing the intensity ratio of the experimental 
and theoretical curves, like in the one-gap approximation. 
This expression was used to fit experimental curves and 
to derive the parameters A 1^2, ri.2 (or 71.2), Z ^ S and K . 
The calculation technique is detailed in j2J|, Appendix. 
The average gap is found from the expression 



AiiC-h A2(l - is:) 



GBTK model 



(3) 



The temperature dependencies of the larger, smaller 
and the average gaps obtained in the two-gap approx- 
imation within the GBTK model with a fixed contri- 
bution K are shown in Fig.[TT] (also see Appendix B). 
Although K is fixed, the critical temperatures obtained 



through a BCS extrapolation are difi'erent for the larger 
and smaller gaps in the PM region, which may point to 
the two-gap character of superconductivity with a weak 
interband scattering. Such behavior is impossible in the 
case of ordinary OP anisotropy where these is only one 
critical temperature, according to Pokrovsky's theorem 
Besides, BCS extrapolation can be employed to 
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estimate Tc in the AFM region for the larger gap in the 
a6-plane and for the smaller one in the c-direction. The 
temperature dependencies of the broadening parameters 
r are illustrated in Fig.IH In the AFM region Pi > P2 
in the a6-plane and Pi P2 in the c-direction. The 
comparison with the one-gap calculation (Fig.[7l) shows 
that the shapes of the curves in Fig. [7] are closer to P2 , 
i.e. the smearing of the high-energy part of the gap is 
important in the one-gap approximation. The temper- 
ature dependencies of gaps and broadening parameters 
discussed in this section are of illustrative character be- 
cause of i4r=const. The goal was to show that they have 
different shapes and Tc differs from BCS-extrapolated 
critical temperatures. Note that the average gaps shown 
in this figure practically coincide with those calculated in 
the one-gap approximation (Fig.[S]). 
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FIG. 11: Temperature dependencies of gaps A (GBTK [H 
[21]) calculated in the two-gap approximation by Eq.([2ll for 
the PC in Fig.[T] The average gaps are Aavor = AiJC'-f A2(l — 
K). The polynomial fit runs through the points outside the 
BCS approximations. The scaling factors are Sab~0.35 and 
Sc=0.5. The contribution to conductivity from the smaller 
gap is K—0.8 (if=const) 



BB model 

The two-band approximation in the BB model was 
considered in [l^ using a free scaling factor S. Here 
we report the results obtained with a fixed scaling fac- 
tor (S'=0.25) for both directions (see Appendix B and 
Fig.BU]). In contrast to the GBTK approximation, the 
contribution K of the smaller gap to conductivity was 
not constant. Nevertheless the number of fitting param- 
eters was the same because at T >2 K 71 turns zero in 
both directions and hence is no longer a fitting param- 
eter. Besides, unlike the broadening parameter F, the 
pair- breaking parameter 7 does not increase the uncer- 
tainty in calculating the r.m.s. deviation (cf. Fig.fTSlwith 
free S and Fig. [191 Appendix). 

The temperature dependencies of the larger A2, 
smaller Ai and the average OP Aavcr, Eq.Q, are shown 
in Fig.[ll 

Note that the use of another model and a non-fixed 
K to conductivity affect the behavior of the temperature 
dependencies of the OPs (mainly the smaller OP) in com- 



FIG. 12: Temperature dependencies of broadening parame- 
ters r calculated in the two-gap GBTK model. r2 and Fi 
are for the larger (A2) and smaller (Ai) gaps, respectively 
fFig.llip. For clarity, a polynomial fit is drawn through the 
points 



parison to the GBTK model and BB approximation with 
varying S [l^. In the a&-plane the smaller OP decreases 
rapidly with temperature and the BCS-extrapolation 
gives Tc -3.3 K. At T >3 K Ai in the AFM region 
changes into the BCS-dependence with Tc ~6.15 K. In 
the PM region there is a region of a smoothly decreas- 
ing gap which persists up to the normal state. This re- 
gion may be due to the interband interaction. In the 
c-direction the BCS-extrapolated Tc values for Ai are 
somewhat higher: -4.8 K (AFM region) and -8 K (PM 
region) . 

The temperature dependencies of the larger OP A2 es- 
pecially in the PM region closely resemble the shapes of 
the curves for the OP calculated in the one-gap approxi- 
mation within the BB model with a varying scaling factor 
(Fig.|9]). The reason for the coincidence is quite obvious. 
In the PM region the shape of the curve is mainly depen- 
dent on the larger OP because its value is several times 
higher than that of the smaller OP. In contrast to the 
one-gap approximation, the contribution of the smaller 
OP to the conductivity holds the scaling factor S con- 
stant. Since in the low temperature region OPs of differ- 
ent bands are comparable in magnitude, the distinctions 
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FIG. 13: Temperature dependencies of OPs A (BB model 
[2^) calculated in the two-gap approximation by Eq.© for 
the PC in Fig.[l] The polynomial fit is drawn through points 
outside the BCS approximation. 5=0.25 in both cases 

between the larger OP of two-band calculations and the 
OP in the one-gap approximation are more explicit in the 
shapes of the curves. The critical temperatures obtained 
by BCS-extrapolation for A2 in the PM region practi- 
cally coincide with the one-gap calculation (Tc ~12 K) 
and exceed the superconducting transition temperature 
of the compound. On approaching T ~ 10 K the OPs 
start to depart from the BCS- dependence tending to 
at Tc of the sample. It is interesting that the tempera- 
ture dependencies of the average OP Aavcr have sections 
in both directions in the PM region that decrease almost 
linearly and go through zero at Tc of the sample. 

Let us consider the temperature-dependent contribu- 
tion K to conductivity made by the smaller gap (Fig.fTH 
and Eq.©). 

At low temperatures the dependence K{T) is readily 
predictable qualitatively: K decreases with a decrease 
in Ai. On a further rise of temperature these param- 
eters exhibit a correlated change in the c-direction: the 
growth of Ai due to the reduced influence of spin density 
waves is attended with an increase in K . The somewhat 
different degrees of the changes in these parameters in 
the c-direction and the absence of a similar correlation in 
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FIG. 14: Temperature dependencies of the contribution to 
conductivity from the smaller OP fFig. 1151 1 calculated in the 
two-gap BB model 



the a6-plane might be attributed to the increasing pair- 
breaking parameter for A2 (Fig. ll5|) but this assumption 
is in conflict with the sharp growth of K in PM region. 
It is then, reasonable to assume that K grows because 
the relative share of the FS containing the large gap de- 
creases. This fact can be explained as follows. 

It is pointed above that the scaling factor S is depen- 
dent in particular on what part of the whole FS is oc- 
cupied by superconductivity. If S'=const., this share is 
independent of temperature. Also, the areas of the FS 
bands that contain a larger and a smaller gap are also 
temperature-independent. This suggests that if the FS 
area with a large OP decreases, superconductivity in this 
part of the band is not suppressed fully. In terms of our 
assumption, the OP in the "vacant" part of the band re- 
duces and becomes comparable to the OP in the second 
zone. We may thus conclude that the superconducting 
share of the total FS area remains invariable. This can 
account for the redistribution of the relative FS shares 
between the large and small OPs. Near Tc on approach 
oiT ^ IQ K the contribution of the large OP to conduc- 
tivity falls below ~20%. As a result, A2 deviates from the 
BCS-dependence and turns rapidly to zero. The physi- 
cal reason for the reduction of the FS area with a large 
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FIG. 15: Temperature dependencies of the pair-breaking pa- 
rameter 7 calculated in the two-gap BB model. 72 and 71 
correspond to the larger OP (A2) and to the smaller OP (Ai) 
respectively (Fig. ll3|) . For clarity, a linear fit is drawn through 
the points 



OP may be connected with the spin fluctuations which 
enhance with temperature. It is likely that larger OP oc- 
curs in the FS part corresponding to the region in which 
the crystal lattice has magnetic moments. This assump- 
tion is supported by the values of pair-breaking param- 
eters 72 of large OP which increase with temperature. 
Note that the 72 magnitudes of larger OP calculated in 
the two-gap approximation within the BB model for non- 
magnetic LuNi2B2C are lower and decrease (faster in the 
a6- plane) with temperature (see [2^, Fig. 9). Since spin- 
flip scattering is absent in this compound, the parameter 
7 in this calculation tends to exceed the degree of the 
superconducting gap broadening. Note that in YNi2B2C 
the gap is most broadened in the low-energy part {29], 
Fig. 6). Assuming that this reasoning holds for ErNi2B2C 
too, we can conclude that no broadening of larger and 
smaller gaps occurs above 2 K (Fig.[T5]) where the param- 
eter 7 is determined solely by the processes of scattering 
at magnetic moments. It is therefore most justified to 
apply the two-gap modification of the BB model in this 
case. The temperature dependencies of broadening pa- 
rameters 7 are shown in Fig. 1151 Two practically parallel 
parts of a linear growth of 72 in the AFM and PM re- 



gions are distinctly seen in the c-direction. There is only 
one linear portion in the a6-plane, which may indicate 
that spin-density waves are ineffective during a magnetic 
transition. Note that the illustrated curves correlate to 
a certain degree with the curves describing the contribu- 
tion of a smaller OP to conductivity (Fig.lTl)). 




T, K 

FIG. 16: The temperature dependence of the lager gap A02 ~ 
A2 (1 - 72''^) (BB model) for the ErNi2B2C-Ag PC in the 

a6-plane and the c-direction (Fig.ll3[). Solid lines are BCS 
extrapolations. Since 71 =0, the smaller gap behaves similarly 
to the small OP. 



Fig. mi illustrates the temperature dependence of the 
energy gap Aq corresponding to A2 of the larger OP. 
There are different BCS extrapolations for changing to 
gapless superconductivity in the AFM (9 K and ~6 K 
in the a6-plane and c-direction, respectively) and PM 
(^11 K in both directions) regions. The gap, like the 
OP, grows during the transition to the PM state in the 
c-direction and has a monotonic dependence in the ab- 
plane. 

It is interesting that the tunnel parameter Z (Fig. fTT)) 
has a feature at Tn only in the c-direction, like in the 
one-gap case (Fig. [5]), which may be attributed to the 
effect of spin-density waves. 
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FIG. 17: Temperature dependencies of the tunnel parameter 
Z (two-gap BB approximation for the ErNi2B2C-Ag PC in 
the a&-plane and the c-direction) . 



CONCLUSIONS 

A detailed analysis of the temperature dependencies of 
PC Andreev reflection spectra dV/dI{V) has been per- 
formed for ErNi2B2C (Tc «11K) in the a6-plane and 
c-direction using one-gap and two-gap approximations. 
Two models were used: the traditional GBTK model 
including the broadening parameter F 21| and the BB 
model [20| in which the parameter 7 characterizes the 



pair-breaking effect of magnetic moments (likely Er). For 
the first time the calculation has been made comparing 
both the shape and the intensity of experimental and the- 
oretical curves. This has decreased the degree of uncer- 
tainty in the temperature dependence A(T) for contacts 
with high broadening parameters F. 

The following conclusions have been drawn. 

1. An anisotropic effect of AFM ordering has been 
detected irrespective of the data processing model. 
For example, the magnitude of the superconducting 
gap calculated within the one-gap GBTK model 
decreases on transition to the AFM state by ~ 25% 
in the c-direction and by ~ 4-7% in the a6-plane, 
which correlates with the behavior of the averaged 
gap in the two-gap approximation within the BB 



model [19|. 

2. The intensity of the PC spectra dV/dI{V) changes 
in correlation with the gap: the intensity decreases 
monotonically in the a6-plane and has an extremum 
near the AFM transition in the c-direction. This 
behavior may be due to the orientation-dependent 
pair-breaking effect of spin-density waves. Thus 
it has been found unambiguously that the AFM 
transition has an anisotropic effect on the super- 
conducting state. 

3. As in [ig|, the pair-breaking parameter 7 increases 
in the vicinity of magnetic transitions, which is nat- 
ural to attribute to the effect of spin fluctuations 
under a change of the magnetic order. 

4. It has been shown that the proper choice of the 
scaling factor S in the one-gap GBTK calcula- 
tion gives the ratio 2A(0)/kTc ^^3.52 for the gap 
in the PM region and its BCS-like temperature 
dependence. The ratio obtained in the one-gap 
calculation within the BB model is 2A(0)/kTc ^ 
4.08 -=-5.14. 

5. The analysis of the models shows that the two-gap 
calculation in the BB model with a fixed scaling 
factor S provides the most adequate information. 
This calculation gives different BCS-extrapolation 
data for Tc of the larger and smaller OPs, which 
points to the multiband nature of superconductiv- 
ity in ErNi2B2C. Its physical sense is that there are 
FS parts with a weaker electron-phonon interaction 
(EPI) in which the temperature induced suppres- 
sion of superconductivity is faster. This is sup- 
ported by the calculation of the anisotropy of the 
EPI parameter in LuNi2B2C which can vary 
from 0.3-^0.8 on a spheroidal FS to 1.0-^2.7 on a 
cube-like FS. 
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APPENDIX. 
Some details of calculation technique 

The technique used (detailed in (i^l) is based on a se- 
lection of parameters that can ensure the smallest r.m.s. 
divergence between experimental and theoretical curves. 
First, the curves dV/dI(V) were normalized to the curve 
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dV/dI{V) taken above Tc and symmetrized. The curve 
fitting was performed in the interval ±8 mV to avoid the 
effect of the phonon features (inflection points) in the 
vicinity of 10 mV in some curves (see Fig.[T]). 



Choice of proper scaling factor 5* 

Theoretically S—l. But S <1 happens very often too. 
The reasons may be as follows. 

1. The electron transport through a PC deviates from 
the ballistic conditions [39l |. 

2. An inhomogeneity at which the PC region at the 
S-electrode side is not fully superconducting be- 
cause of the normal region (regions) near the N — S 
boundary in the superconductor. 

3. A superconducting gap can occur only in a part of 
the FS. For example, the AFM molecular field in- 
duces a gap in some FS parts in ErNi2B2C and sup- 
presses the superconducting gap in these regions. 




FIG. 18: The dependence of the r.m.s. deviation F 
of the shape of the theoretical curve from experimen- 
tal A at 7K in the c-direction. The intensity ra- 
tio S={dV/dI{V))e^p/{dV/dI{V))theoT between experimental 
and theoretical curves is changing from 3.58 to 0.23. The 
lowest F corresponds to the best coincidence of the exper- 
imental and theoretical curves is obtained at A=0.96 meV 
and 5=1.3. The narrow parabolic error curves were taken at 
S^const (one-gap GBTK model [IH 

A more exotic case of >1 is also possible. For exam- 
ple, to improve the description of the shapes of smeared 
experimental curves, the parameter F which assumes a 
finite lifetime of carriers is usually increased. As a result, 
the intensity of theoretical curves decreases. However, if 
the spectra dV/dI{V) are smeared e.g., because of dif- 
ferent values (distribution) of the superconducting gap 



in different FS parts due to anisotropy or multiband su- 
perconductivity, the theoretical F-broadening curve can 
approach the shape of the experimental curve, but its in- 
tensity will be lower, and the scaling factor S can exceed 
1. This can be used as a criterion of validity for a model. 
As an example, we consider the spectrum from Fig. 14 in 
[23 ]. The gap is distributed in the interval 1-3.35 meV 
(insert). The calculation in [1^ gives: 

1. GBTK (1-gap): A=2.565 meV, F=0.523 meV, 
Z=0.8, 5=1.467. 

2. GBTK (2-gaps): Ai=2.18 meV; A2=2.99 meV; 
Fi=0.36 meV; F2=0.067 meV; Z=0.78; K=0.56; 
Aavor=2.538 meV; S'=1.144. 



3. BB (1-gap) 
S'=0.968. 



A=2.79 meV, 7=0.04, Z=0.74, 



4. BB (2-gaps): Ai=2.15 meV; A2=2.984 meV; 
71=0.046; 72=0.011; Z=0.75; X=0.315; 
Aavcr==2.72 meV; S'=0.968. 

Thus, the one-gap GBTK calculation results in the high- 
est value in the 5' estimate. Therefore one-gap GBTK 
calculations [i^l can be regarded as oversimplified and 
can be used only as a first approximation. In ordinary 



superconductors, e.g., Zn [4l|, S ^1 is independent of 



temperature because the gap opened isotropically on the 
whole FS at T <Tc. 

In Fig. 18 the smooth broad arc-like curve corresponds 
to the best coincidence of the shapes of theoretical 
and experimental curves for the gap A in the interval 
0.6 2 meV. The lowest error is at A=0.96meV. The 
scaling factor varies along the curve from S =3.58 at 
A=0.6meV to S'=0.23 at A=2meV and S = 1.3 at the 
lowest F. S' >1 is possible only assuming the gap dis- 
tribution (see above), and therefore the lowest F (shape 
error) alone is not sufficient to be a criterion. Besides, 
for the curves taken at higher F values a comparatively 
small change in the shape of the temperature- neighboring 
curves dV/dI{V) can shift arbitrarily the error minimum 
and the corresponding A. 

Unlike the GBTK model, the lowest-error curves ob- 
tained in the BB model have distinct minima (Fig. [T9| . 
In our case the criterion of the proper choice of the scaling 
factor S was its value at which the shapes of the theo- 
retical and experimental curves coincided most closely in 
the whole interval of temperatures used. To avoid over- 
loading. Fig. 20 contains the F and S calculation for two 
characteristic temperatures: T=3.5 K in the middle of 
the AFM region and T—7 K, i.e. above the tempera- 
ture of AFM ordering. It follows from Fig. 20 that S can 
vary from 1 to 0.2, being the lowest at F — >-0 fFig.lSO]). 
The vertical lines show the best 5'-values corresponding 
to the minimum error along with their associated gaps 
(OPs) and broadening (pair-breaking) parameters. Note 
that in the BB model a departure from the best-chosen 
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FIG. 19: The dependence of the mutual r.m.s. deviation F 
of the theoretical and experimental curve shapes from A (BB 
model) in the c-direction at different temperatures. S varies 
along the curves and at the minima, T=3.5, 5.5 and 7 K, 
being 5=0.24, 0.2 and 0.22, respectively 



5'- value caused a sharp increase in the error F, and there- 
fore the limits of S variation are rather narrow. On the 



other hand, the GBTK model [21| allows more freedom 
of selecting S, especially in the c-direction. However, this 
variation has practically no effect on the shape of tem- 
perature dependence A(T) because the curves A(S') are 
nearly parallel at both temperatures. The gap can be es- 
timated readily at another S value using the dependence 
A(5). 



One-gap BB calculation with S^const 

Unlike the GBTK model, the lowest-error curves ob- 
tained in the BB model have distinct minima (Fig. [191). 
Using a fixed S impairs considerably the quality of fitting 
in this approximation. The results of the one-gap calcu- 
lation allowing for the minima in the curves are illus- 
trated in Fig. [9] and the dependencies S{T) are shown in 
Fig. 131 The decrease in the S'-value is due to the two-gap 
character of superconductivity in this compound, which 
determines different OP magnitudes and a different (in 
comparison to the GBTK result) shape of the density of 
states. The relative share of the FS with a larger OP 
decreases at rising temperature (see the Two-Gap BB 
calculation section). Used in BB theory the density of 
states terminates abruptly at the gap edge. In the one- 
gap approximation of the experimental curve this entails 
termination of the contribution to conductivity from the 
FS parts with smaller OPs. However, since the spectrum 
intensity is estimated assuming that the OP obtained is 
related to the whole Fermi surface, S decreases. The dif- 
ference between the dependencies S{T) in the c-direction 
and the a6-plane can be attribute to the anisotropic ef- 



fect of spin-density waves (see the above interpretation of 
the quantity M and Fig.jSj. The temperature-dependent 
AFM molecular field induces a gap in some parts of FS 
in the c-direction and destroys the superconducting gap 
in them. An AFM gap is more probable in the FS parts 
which have a smaller superconducting gap taking into 
account the temperature dependence of the OP in the 
corresponding temperature region. 

The sharp decrease in the OP on approaching Tc can 
also be caused by the two-band structure of the FS and 
the reduction of the FS share with a larger OP. There 
may exist a certain minimal FS share with a large OP 
below which superconductivity is destroyed rapidly due 
to the interband interaction. In our case the departure 
from the BCS-dependence started in both directions at 
T ~10K at which the contribution of the larger OP to 
conductivity dropped below 20% (two-gap calculation). 



Reducing the number of fitting parameters for the 
two-gap GBTK model 



As noted in the Introduction, in the general case the 
two- gap calculation by Eq.([2l) involves seven fitting pa- 
rameters (Ai, A2, Fi, Z, S and K). This has an 
unfavorable effect on the estimates obtained. It is there- 
fore desirable to minimize the number of such parameters 
on the basis of physically reasonable limitations. In the 
strict sense, the scaling factor S is not a fitting parameter. 
It does not enter into the theoretical formulas describing 
current-voltage characteristics. It is intended to equalize 
the intensities of experimental and theoretical curves on 
calculating the mutual r.m.s. deviation of their shapes. 
However, a use of S= const restricts the range of permis- 
sible values for the rest of the fitting parameters. A fixed 
S reduces the number of fitting parameters at least by 
one parameter fFig.lTSl). We did not fix the tunnel pa- 
rameter Z which was found to be approximately constant 
except near Tc. As a result, we have four fitting param- 
eters instead of seven. Besides, it is reasonable to fix the 
relative contribution K to conductivity from each band. 

Thus, the calculation was made using fixed 5=0.35 
in the a&-plane and 5=0.5 in the c-direction. The rela- 
tive contribution of the smaller gap to conductivity was 
K=Q.S in both directions. Fixed K actually couples A 
and F, which can distort their temperature dependencies. 
Indeed, if, for example, the first gap remains constant and 
the second one decreases, the relative contribution of the 
first gap increases (it is assumed that the FS share of 
each gap does not change). To exclude this, it is neces- 
sary to decrease the broadening parameter of the second 
gap or to increase it for the first one. Identical temper- 
ature dependencies for gaps and broadening parameters 
is the simplest version free of distortions. 
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FIG. 20: Dependence of the relative deviation of the shapes F of theoretical and experimental curves from S calculated in the 
GBTK (a, b) and BB (c, d) models at r=3.5 K and 7 K in the a6-plane and in the c-direction. The lower part of each figure 
shows the corresponding dependencies of the gap (OP) A and the broadening F (pair-breaking 7) parameters. The vertical 
line marks the best choice of S corresponding to the closest coincidence of the curve shapes at both temperatures 



NOTE ADDED IN PROOF 



tribute to DoS and therefore our two-band approach is 
reasonable. 



After the paper was sent to arXiv, we were aware 
about recent three-dimensional study of the Fermi sur- 
face of LuNi2B2C in This study shows that I) the 
Fermi surface topology of the rare-earth nickel borocar- 
bides varies little for rare-earth elements such as Er, Tm 
and Yb, 2) there are 3 bands which contribution to the 
density-of-states (DoS) at the Fermi energy is 0.24%, 
22.64% and 77.1%. That is, two bands basically con- 
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